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O '. Abstract 

Time analysis of oscillations of a particle between wells in the one-dimensional 
double-well potential with infinite high outside walls, based on wave packet use 
and energy spectrum analysis, is presented. For the double-well potential of the 
form x 2 + 1/x 2 in the external regions, an exact analytical solution of the energy 
spectrum is found (by standard QM approach), an analysis of oscillation periodicity 
is fulfilled, an approach for exact analitical calculation of the oscillation period is 
proposed (for the first time). 
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1 Introduction 



A quantum system, representing a particle in the potential field in the form 
of two wells with a finite high barrier and infinite high outside walls, is used 
in many different tasks of physics and chemistry [1]. 

If the tunneling of the particle through the barrier in such system is possible 
then there are transitions of this particle between the wells, named as "oscilla- 
tions". In the general case, the oscillations are not periodic (and also harmonic) 
motion. For their description one can introduce the following characteristics: 
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• the time duration, after which the system returns into its initial state (this 
time characteristic can be named as oscillation period of the particle between 
the wells [2]) or period of Poincare's cycle); 

• the time duration, after which the system has passed into a state closed as 
much as possible with initial one. 

Instanton methods [3,4] are powerful tools for realization of a time analysis 
of a behaviour of the particle in the double-well potential. An alternative 
and less widespread approach for the time analysis of the double-well system 
behaviour consists in use of wave packets in spite of the fact that such system 
has discrete energy spectrum only [5] . If to localize the packet in the initial time 
moment inside one well than one can suppose that a maximum of this packet 
has passed into another well after some time duration. One can define the 
oscillation period of the particle between wells for the systems with a periodic 
motion and also one can define the time duration of the most probable return 
into the initial state for the system, which motion is not periodic, on the 
basis of transition of the packet maximum or the center of mass of this packet 
between the wells. Note that the packet maximum and the center of mass of 
the packet describe the tunneling process of the packet through the barrier in 
different ways. Here, if in the tunneling the localization of the wave packet in 
the space region of the barrier is improbable (apparently, it is possible only 
in one case of essential contribution into the total packet of its component for 
sub-barrier energies), than a motion of the center of mass of the packet in the 
barrier region looks natural enough. 

An approach for calculation of these time characteristics for the double-well 
systems on the basis of analysis of energy spectrum is presented in this paper. 



2 An analysis of the time periodicity of wave function on the basis 
of the energy spectrum 

Let's consider a system with the discrete energy spectrum where energy levels 
are located one from another at such distances, for which one can calculate 
exactly the largest general divisor: 

E n = E + A • n, (1) 

where n G 0, N (N is a set of natural numbers). In particular, a harmonic 
oscillator, a particle inside box, a charged particle in the Coulomb field satisfy 
to this condition. Wave function of such system has a form [6]: 

*(x,t) = e- lE ^ h J29nMx)e- l(En ~ Eo)t/h = e~ lEot / h J2 9nM^)e~ lAnt/ \ (2) 

n n 
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where ij) n {x) are the orthonormal functions of the stationary states of the 

system satisfied to equation H 4> n {x) = E n ^jj n (x); H is Hamiltonian of the 
system; I] |g n | 2 = 1. We choose the time moment t — for time zero. 

n 

Consider the function 

f(x,t)=9(x,t)e iE ° t ' h (3) 



in the time interval [— ^]. It satisfy to Dirichle's condition [7]: a) one can 
divide this interval into a finite number of intervals where the function f(x, t) 
is continuous, monotonous and finite; b) if to is a point of discontinuity of 
the function f(x,t), then f(x,t + 0) and f(x,t — 0) exist. And the sum 
J2 n 9ni J n(x)e~^ En ~ E °^ t ^ h in (2) is an expansion of the function f(x,t) into the 
Fourier trigonometric series (with constant coefficients A/K) by variable t. 
Such series is convergent at any point t of the interval [—7^,7^], where the 
function f(x,t) is continuous in t. In any break point to the Fourier series 
converges to /( 3: ' f o+ )+/( :E ' f o~ ) _ Therefore, the function f(x,t) is periodic by 
variable t in the interval considered above. This periodic dependence can be 
extended into whole range of definition of this function by t. One can calculate 
the period by such a way: 

T f = ~1T- ( 4 ) 



An exact analytical dependence of the periodic function f(x,t) on time vari- 
able t is determined by the total set of eigenf unctions tp n (x) with coefficients 
g n at any point x and is changed in dependence on x. 

We find out an influence of the ground state with the level E on a period- 
icity of the function ty(x, t) in time. Consider the following example. Let the 
function f(x ,t) at point x = x be harmonic by variable t (assume that this 
function can be written as const ■ e~ lAt ^ h ) : then we obtain: 

V(x , t) = const • e -iE t/n e -iAt/h = congt . e -i(E +A)t/h_ (5) 



One can see that the function ^/(x ,t) is harmonic and periodic by time vari- 
able t also. We find an exact analytical solution for a period of the function 
ty(x 0: t) with taking into account the first level E : 

_ 2nh _ 2nh 



Using (6), we can analyse a contribution of the first level E in the period TV 
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In particular, if to use the parabolic function U(x) = mw 2 x as the potential 
energy of the system (where m is a mass of the particle, w is an oscillation 
frequency of the particle in classical mechanics), then the energy spectrum of 
the system has a form: 

E n = (n + \)hw,n = 0,1,2... (7) 



If there is a coordinate x where the function f(x ,t) is harmonic by time 
variable t then the oscillation period relatively this coordinate decreases in 1, 5 
times with taking into account the first level E = ^ (i. e. the contribution 
of the first level is not small): 

_ 2irh _ 2nh 2nh 2nh _ 2 2irh _ 2 

~ "a" ~ T^T' new ~ e^Ta ~ {\ + i)hw - 3 ' 7^7 ~ 3 old - (8) 



Therefore, the function ^{x, t) is harmonic in time at point x only in such 
case when at this point x the function f(x,t) is harmonic in time also. 

If the function f(x,t) is not harmonic then one can consider the function 
ty(x,t) as periodic in time approximately and evaluate its period. Let fi(t) 
and f 2 (t) be two periodic functions with periods 7\ and T 2 , accordingly. We 
assume that these functions are harmonic in enough small neighborhood of 
the point (x, t). The period T (i. e. the time duration, after which the system 
has passed into a state closed as much as possible with initial one) can be 
calculated from the following relation: 

111 

f ~ 7\ + TV ^ 



Taking into account that the functions f(x,t) and exp (—iE^t/h) are periodic 
by variable t, one can determine conditions of the periodicity of the function 
V(x,t). 

Any quantum system with discrete energy spectrum performs a finite motion. 
But if there is the largest divisor A determined exactly for energy level of 
this system and the condition (1) is satisfied, then in general case one can 
consider the time behaviour of this system as periodic approximately and one 
can calculate the period more accurately by (6) or (9) (i. e. with taking into 
account of the first level) whereas the expression (4) (i. e. without taking into 
account the first level) is less accurate. The calculation accuracy of oscillation 
period can be estimated on the basis of the set of orthogonal functions ip n (x) 
with coefficients g n at selected point x. This method of the estimation of 
the periodicity of the wave function can be used for the time analysis of the 
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oscillations of the particle between wells in the double-well potential with 
infinite high outside walls. 



3 Exactly solvable models with symmetric double-well potential 



Apparently, the spectrum of the double-well potential with any form is not 
equidistant. Nevertheless, we find the double-well potential form, for which 
there is the exact analytical dependence of the wave function on time and one 
can analyse its periodicity on the basis of the energy spectrum. 

Let's consider a quantum system, representing a particle in the symmetric 
double-well potential field of a form: 



U(x) 



raw 



B 2 

(x + Xq) 2 + - — , for x > a > 0; 

(x + XqY 



C - Dx 2 , 

(X - Xq) 2 + 



B 2 



(x - Xq) 



for — a < x < a; 
-, for x < —a < 0; 



(10) 



where m > 0, w > 0, B > 0, C > and D > 0. C determines the barrier 
height. We assume that B, C and D have such values that U(x) is continuous 
at points x = ±a. Minimums of two wells are located at points x\ = —\/B + xo 
and x 2 = \f~B — xq. 

The potential (10) is symmetric. In result of sub-barrier tunneling and above- 
barrier propagation there are transitions of the particle between wells, i.e. its 
oscillations. One can analyse the period of such oscillations on the basis of the 
energy spectrum. For finding the spectrum it is need to resolve the stationary 
Schrodinger equation and take into account all conditions, to which the wave 
function satisfies. In the calculation of the energy spectrum for the potential 
(10) the boundary conditions at x — > — oo and x — > +oo play the important 
role. According with analysis, the conditions of the continuity of the wave 
function and its derivative at points x = ±a do not differ the localization 
of the energy levels. Therefore, the solution of the stationary Schrodinger 
equation in the external regions x < —a and x > a has caused the main 
interest in calculation of the energy spectrum. 
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3.1 Calculations of the energy spectrum 



We use new parameters: 



2 2 

m w 



m 2 w 2 B 2 
W 



h 2 



,K = 



(11) 



Then the stationary Schrodinger equation in the external regions transforms 
into a form: 



d 2 j) 
dx 2 



+ (g + Fx 2 + i) = 0, 



(12) 



where 



x = x — xq, for x < —a; 
x = x + xq, for x > a. 



(13) 



We find a solution of this equation. Fulfill the replacements: 

e = ocx 2 , a = ™ = v=f, m = (I ) 



(14) 



Then the equation (12) transforms into the standard Whittaker form 

d 2 w 



de 



+ 



1 G 1 / 3 K\ 1 



4 4 V /Z F£ V16 4/£ 



w = 0. 



(15) 



Include the following parameters 

G.IK 1 , , ft 

k = -£j=,» =-- T , a =--k + »,c = l + 2» 



(16) 



and fulfill the following replacement 

y(0 = r c/ v/vo- 



(17) 



Then the equation (15) transforms into the hypergeometric equation of the 
form: 



^ 2 + (c-O^-ay = 0. 



(18) 
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Its partial solutions can be presented using the hypergeometric function F(a, c; £) 
in a form: 



yi(0 = F(a,c;S), 

2/ 2 (0 = £ 1 - c F(a-c+l,2-c;0, 

y 3 (S) = e£F(c-a,c; -£), 

y 4 (0 = e 1_c ^(l-a,2-c; -£)■ 



(19) 



Consider the case: c ^ Z. Here, the solutions y 3 and y 4 can be written through 
yi and y2 by Kummer's transformation [8]. Therefore, they depend linearly on 
yi and y2- We obtain the first two solutions y\ and y 2 with initial variables: 



In accordance with the finiteness condition of the wave function, one can write: 



for tp 2 (x) : -a + 2/i G 0, TV; 2/i £ N. 

After satisfying of these conditions the series (which define the hypergeometric 
functions for the solutions ipi and ^2) transform into polynomial, and the 
energy spectrum becomes discrete. One can find from the expressions (21) 
that the solutions tpi and ip 2 can not be used at the same time. But both ip 1 
and ip2 corresponds to the same energy spectrum. Write expressions for the 
energy spectrum and wave function (we point out the wave function in the 
external regions only): 



< ip hn (x) = a^i-^x^-^e-^^Fi-n, 1 - 2/i; ax 2 ), for x<-a andi> a; ( 22 ) 
^2,n{ x ) — a 1 / 2+ ' I x~ 1 e~ ax2 / 2 F(— n, 1 — 2/i; ax 2 ), for x < —a and x > a; 



< ip hn ( x ) = a l > 2+ ^x~ 1+2 ^e- ar£2 l 2 F(-n, 1 + 2/i; ax 2 ), for x<-aandi> a; (23) 



^(x) = ^1/2+^-1+2^-^72^ ( a> 1 + 2/x; ax 2 ) , 
ijj 2 (x) = a 1 / 2+ ^x- l e- as2 l 2 F (a - 2/i, 1 - 2/i; ax 2 ) . 



(20) 



for Vi(z) : a G 0, — iV; 2/i £ —N, 



(21) 





^2,n(^) = a 1/2 ~^x _1 e~' 



l 2 F(-n, 1 + 2/i; ax 2 ), for x < 



a and x > a. 
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Here 



1 / 4m 2 w 2 B 2 | \x - x \, for x < -a; 



\x + x \, for x > a, 



n E 0, iV (N is the natural number set) and the total set of the levels E n 
consists from sets E+ and E~ . 

One can consider a co-existence of two eigenfunctions ipi and ip2 (not depen- 
dent linearly one from another) for any level E n as the doubly degeneracy of 
the energy spectrum. In similar cases, an additional quantum number is used 
for marking the difference between such states. Inclusion of the asymmetry in 
such double-well potential gives the double splitting of the energy spectrum 
and leads to degeneracy removal. Instanton methods give a similar result. 



3.2 The time analysis of the particle oscillations between the wells 



We write the expressions (22) and (23) for the energy spectrum by such a way: 



E^ = 2hw Q ± /ij + 2hwn = Eq + A ■ n, 



E± = 2hw(l±fi), (25) 



2 



A = 2hw. 



From here one can conclude that two sets of the energy levels and E~ 
represent independently equidistant spectrum. One can study them separately 
as described two independent waves with their periods. Without of taking into 
account the first levels E£ and Eq , the periods for these waves are equal: 

_ 2nh 7r , , 

T > = -5" = » (26) 



With taking into account the first levels Eq and E , a dependence of the wave 
function ^(x, t) on time is defined by such a way: 

V(x,t) = e- iE ° t/h f + (x,t) + e- iE « t/h f-(x,t), (27) 



where f + (x,t) and f~(x,t) are the wave functions of these waves without of 
consideration of the levels Eq and Eq . The expression (27) represents the 
exact analytical dependence of the wave function on time. One can consider 
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approximately these two items as periodic with the periods T + and T calcu- 
lated on the basis of (6): 

_ 2nh 2nh _ 2nh 2nh 

~WTa~W ~W+*~W (28) 



One can calculate approximately the period of the particle oscillation between 
the wells in the potential (10): 

111 

- = + . (29) 

T T+ + T~ T+ — T~ 



In accordance with (25), a distance between two nearest levels E+ and E~ 
(considered in some tasks as a splitting value of energy spectrum E n in result 
of tunneling through the barrier) can be calculated by such a way: 

AE = El - E- = Hw^l + ^B 2 . (30) 



If the barrier form satisfies to conditions of use of semi-classical methods then 
one can find a dependence of the penetrability coefficient D of the barrier on 
the oscillation period T and on the largest divisor A [2]: 



4m 2 w 2 „ 9 
£>~7T 2 (l + — —B 2 

7r 2 16m 2 B 2 



7T 2 (1 



m 2 B< 



A 2 = 



7T 2 1 + 



h 2 



J^2 



(31) 



4 Conclusions and perspectives 



The approach for the time analysis of the double-well systems on the basis 
of energy spectrum is presented in this paper. For the double-well symmetric 
potential with the form x 2 + B 2 /x 2 in the external regions the energy spectrum 
is calculated exactly analytically and the analysis of the time periodicity is 
fulfilled. The approach to approximate calculation of the oscillation period is 
described. Similar potentials were studied in [9] . 

During last two decades an essential progress in study of quantum system 
properties has achieved after rapid development of methods of supersymme- 
try in their application to quantum mechanics (here, one can note a review 
[10], which should be the best SUSY QM review at opinion of the author). 
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Application of such methods (in particular, see last subsections in [11]) to the 
analysis of periodicity of the particle motion (oscillations) between two wells 
in the double-well potential looks perspective and interesting in future study 
of periodical systems. For example, a question about change of the period- 
icity characteristics of the particle oscillations after going to new isospectral 
potentials from the potential considered in this paper can be interesting. 

For the double-well non-periodic quantum system one can select "quasi-cycles" 
(after which the system has passed into a state closed as much as possible with 
initial one) with needed accuracy [2]. The maximal values of wave function of 
such systems are localized inside a finite space region and one can calculate 
the oscillation period with needed accuracy (with taking into account the 
needed number of "quasi-cycles"). From the other side, it is interesting to 
use such approach in the generalization of the SUSY QM methods, developed 
for the analysis of the periodical quantum systems (i. e. for the analysis of 
quasi-periodical systems). 
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